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Abstract 

We study the Kowalevski expansions near singularities of the swinging At- 
wood's machine. We show that there is a infinite number of mass ratios M/m 
where such expansions exist with the maximal number of arbitrary constants. 
These expansions are of the so-called weak Painleve type. However, in view of 
these expansions, it is not possible to distinguish between integrable and non 
integrable cases. 
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1 Introduction 



The swinging Atwood's machine is a variable length pendulum of mass m on 
the left, and a non swinging mass M on the right, tied together by a string, in 
a constant gravitational field, see Figure ([I]). The coupling of the two masses 
is expressed by the fact that the length of the string is fixed: 
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^/x^ + y^ +\z\ = L, ^ x'^ + y"^ = {\z\ - L) 

Up to a choice of origin for z, one can assume L = 0, so the constraint is 
the cone z'^ = x"^ + y"^. To describe the dynamics we choose to work with 
constrained variables and write a Lagrangian 

where A, a Lagrange multiplier (of dimension MT^"^), has been introduced, 
whose equation of motion enforces the constraint. The equations of motion 
read : 

mx = Xx (1) 
mij = —mg + Ay (2) 
Mz = -Mg - Xz (3) 
= x'^ + y^-z^ (4) 

From these equations one can express A in terms of positions and velocities: 

XX + yij — zz + g{y — z) mM z^ — x^ — ij'^ + g{y — z) 

Alternatively, rescaling 

1 1 1 

X — > —^x, y — > —^y, z z 

Jm Jm JM 



we can view the system as a unit mass particle moving on a cone 

z = — (x^ + y ) 
m 

subjected to a constant field force 



















\-g^j 
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Figure 1: Swinging Atwood's machine. 



The slope of the force in the (y, z) plane coincides with the angle of the cone. 

The swinging Atwood's machine has been studied in great detail by N. 
Tufillaro and his coworkers, see [1-7]. They have first studied numerically the 
equations of motion and shown that for most values of the mass ratio M/m 
the motion appears to be chaotic, however for some values, like 3, 15, etc. 
the motion seems less chaotic and could perhaps be integrable. In a further 
study, Tufillaro [4] showed that the system is indeed integrable for M/m = 3 by 
exhibiting a change of coordinates, somewhat related to parabolic coordinates, 
in which separation of variables occurs. He was then able to solve the equations 
of motion in terms of elliptic functions, which is quite peculiar since in general 
integrable systems with two degrees of freedom can be solved only in terms of 
hyperelliptic functions, such as for the Kowalevski top [8]. He also obtained 
the second conserved quantity which ensures integrability. In the same paper, 
he conjectured that the system is integrable for M/m = 15, • • • ,4n^ — 1, with 
n integer. 

However, later on, Casasayas, Nunes and Tufillaro proved [6] that the sys- 
tem can be integrable for discrete values of the ratio M/m only in the interval 
]1,3], using non integrability theorems developed by Yoshida |9j and Ziglin. 
The essence of the Yoshida-Ziglin argument is to study the monodromy de- 
veloped by Jacobi variations around an exact solution, when the time variable 
describes a loop in the complex plane. The monodromy must preserve con- 
served quantities, but this is impossible in general if the monodromy group 
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is not abelian. In the case at hand one can compute monodromies from hy- 
pergeometric equations and conclude. We have also been informed by private 
communication of J. P. Ramis, that he and his coworkers have proven that the 
swinging Atwood's machine is never integrable except for M/m = 3, using 
methods from differential Galois theory. 

The aim of our paper is to work out the Kowalevski analysis for this model. 
Let us recall the idea of the Kowalevski method. If a dynamical system is 
algebraically integrable one can expect to obtain expressions for the dynamical 
variables in terms of quotients of theta functions defined on the Jacobian of 
some algebraic curve of genus g, where g=2 for a system with 2 degrees of 
freedom. Only quotients may appear because theta functions have monodromy 
on the Jacobian torus, which needs to cancel. Hence denominators which can 
vanish for any given initial conditions and for some finite value, in general 
complex, of time will appear in the solution. Hence the equations of motion 
must admit Laurent solutions -that is divergent for some value of time, with as 
many parameters as there are initial conditions. S. Kowalevski first noted [8], 
that this imposes strong constraints on these equations, from which she was 
able to deduce the celebrated Kowalevski case of the top equation. 

Looking for Laurent solutions to the swinging Atwood's machine equations 
of motion in the integrable case M/m = 3 we first noted that there are none, 
but there exists so-called weak Painleve solutions, that is Laurent developments 
not in the time variable t but in some radical t^l^, generally called Puiseux 
expansions. 

It had already been discovered by A. Ramani and coworkers [1^ that some 
integrable systems require weakening the Kowalevski-Painleve analysis to ob- 
tain expansions at infinity of dynamical variables. This may be explained in 
general, and is certainly the case for our example, by the fact that there is a 
"better" variable which has true Laurent expansions and time itself can be ex- 
pressed in terms of this variable through an algebraic equation which happens 
to produce the given radicals. Moreover Ramani et al. advocated the idea that 
the existence of weak Painleve solutions is a criterion of integr ability, like in 
the Kowalevski's case. 

For our model of the swinging Atwood's machine, we find that there are 
weak Painleve solutions not only when M/m = 15 but for a whole host of other 
values of the mass ratio, all of them corresponding to obviously non integrable 
cases. Hence this model provides a large number of counterexamples to the 
above idea. We then study in detail the solutions around infinity which can be 
extracted from these Kowalevski developments. Using Pade approximants we 
are able to extend these solutions beyond the first new singularity and observe 
how the new singularities obey Kowalevski exponents. 

We also comment on the Poisson structure of the model, which is interest- 
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ing due to the constraints between the dynamical variables, and the Poisson 
brackets of the variables appearing in the Laurent series, which happens to be 
of a nice canonical form. We notice that this illustrates the fact that it is the 
global character of the conserved quantities that is of importance in defining 
an integrable system. 

One of us (M.T.) is happy to acknowledge useful conversations with J. P. 
Ramis and J. Sauloy from Toulouse University, about their work on differential 
Galois theory applied to the swinging Atwood's machine. Finally we are happy 
to thank the Maxima tearG0 for their software, with which we have performed 
the computations in this paper. 



2 Hamiltonian setup. 

The description we have given of the swinging Atwood's machine is a con- 
strained system in the Lagrange formulation, so that the equations of motion 
take a nice algebraic form. 

In the articles [1-7] polar coordinates are used, so the constraint is "solved" 
but the price to pay is the use of trigonometric functions. Using polar coordi- 
nates X = rsin^, y = —rcosO the Hamiltonian reads: 

H = ^^J:^/r + ^,pI +9riM-m cos 9) (6) 

where pr = {m + M)r and po = mr'^9. 

We now give a Hamiltonian description of this system, using as dynamical 
variables the three coordinates x,y,z and the three momenta Px:Py,Pz with 
canonical Poisson brackets. The constraint 

Ci = z^-x'^-y^ = (7) 

generates the flow: 

{CuPcc} = -2x, {Ci,py} = -2y, {Ci,p,} = 2z (8) 

which is also generated by the one parameter group acting on phase space by: 
{x,y,z) {x,y,z), {px,Py,Pz) ^ {Px - l^x,py - fiy,pz + fiz) where /i is the 
group parameter. 

We want to describe the dynamics of our model as a Hamiltonian system 
obtained by reduction of an invariant system under this group action In 

http:/ /maxima. sourceforge.net/ 
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order to do that, consider the functions: 



Ax = zpy + yp^ 

Ay = + XPz 

A-z = xpy - ypx 

These functions Poisson commute with the constraint Ci hence are invariant 
under the group action. They are not independent however, since they are 
related by: 

yAy - xAx + zAz = (9) 
It is easy to check the Poisson brackets: 

{A-XiAyy = —Az, {Ax,Az} = —Ay, {AyjA^} = A^ 

{Ax,x} = 0, {Ax,y} = z, {Ax,z} = y 

{Ay,x} = z, {Ay,y} = 0, {Ay,z} = x 

{Az,x} = -y, {Az,y} = x, {Az,z} = 
Let us consider the invariant Hamiltonian: 
1 



H = 



2(m + M)z2 



AI+AI + -AI 



+ Mgz + mgy (10) 



To check that H generates the equations of motion on the reduced system, 
we compute: 

x = {H,x} = -—jT^izAy yAz] 11 

m + M z"^ \ m J 

y = {H,y} = -——-IzAx + —xAA 12 
m + M z"^ \ m J 

= ^i^^^^'^^'^y^'^ ^^^^ 

The right hand sides of these equations are hnear in the momenta px, Py, Pz, 
however we cannot invert the system uniquely in order to express the momenta 
in terms of the velocities. This is because, due to the symmetry {{H, Ci} = 0) 
we have xx + yy = zz so the equations are not independent. The solution is: 

Px = mx + iix (14) 
Py = my + fiy (15) 
Pz = Mz-jiz (16) 
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where fi is arbitrary. Similarly we compute x = {H,x}, etc... where x, etc... 
are the right hand sides of the above equations. Performing this calculation 
and using the constraint Ci and eq.lIU), we obtain the Lagrangian equation of 
motion ([iHS]), with A given by: 



A 



mM 1 
m + M 



giy 



1 



This coincides with eq.(l5|), as can be cheked using again eqs. (fT4HT6]l and the 
constraint Ci, to express Az in terms of x, y, z. 

Finally we express the energy in terms of velocities still using the con- 
straints. We find: 



E 



m 



which agrees with what we expect from the Lagrangian formulation. 



3 The integrable case. 

In order to understand what sort of Laurent expansions appears in the model 
it is useful to first consider the case M/m = 3 which has been integrated by 
Tufillaro [4l|. Let us recall some of his results. He discovered that using polar 
coordinates {r,9) such that x = rsinS, y = —rcosO and r = z, and setting: 

^2 = z[l + sin(0/2)], ^ _ sin(0/2)] 

then the Hamilton-Jacobi equation separates in the variables rj). These look 
like parabolic coordinates, except that the half-angle 6/2 is used. Knowing ^ 
and rj one can recover x and y by: 

In fact, just for M = 3m, two terms involving couplings between ^ and rj 
disappear, and one gets, with momenta = 4,^(^^ + ff) etc. the expression of 
the Hamiltonian, in which we have set m = 1: 

H = [{pI +pI)/S + 2g{e + il')]/{e + rj') 

Then it is clear that in this case the action S separates as a sum 5'^(^)+5^(r?) 
where 5^ and 5^ obey different elliptic equations (corresponding to different 
elliptic moduli): 

{d^S^f = -Wg^^ + SEf + I^P+iO (19) 
{d^Sr^f = -16grj^ + 8Er]^-I = P^{rj) (20) 
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where I is the separation constant. It can be expressed in terms of dynamical 
variables by substracting the above two equations multiplied resp. by rf' and 
which eliminates E. Moreover we replace: 

We get: 

//16 = ie + v'Kv'f - ee) + a^Wie - v')/{e + v') 

Returning to polar coordinates the integral of motion takes the form: 

I/IQ = r^e\f cos(e/2) - y sin(^/2)] + gr^ sin(^/2) zo^^iQjl) 

We want to see if the equations of motion admit a solution which diverges 
at finite time, and in that case what is the behavior of the Laurent expansion. 
Tiie general solution of the Hamilton-Jacobi equation is: 



dr] 



According to the general theory we get the solution of the equations of motion 
by writing OeS = ce and djS = cj for two constants ce and cj. For c/ 7^ we 
get: 

4£2 m 42 

t + CE = / , di+ / , ' dr] (22) 

" - -J 7m'^--J Tm'" '''' 

For / = 0, the elliptic integrals degenerate to trigonometric ones. We get: 
t + CE = (Vl - + - W) , a = 25/^, uj = g/V2E (24) 



, = K^ ^ V^_^ K2 = ^cr ^25) 

l + 1 + Vl-ar/2' 

so that setting ^ = sin((^^)/-\/a, r) = sin(0^)/-\/a the second equality reads: 

tan((^^/2) = i^tan(0,,/2) 
Using the variable s = tan(0^/2), ^ and ry can be expressed rationally: 

1 2s _ 1 2Ks 
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Finally, one gets the time variation of S = by using eq. (p4|l which implies 



ujdt = dS 



1 



+ 



(1 + 5)2 {K'^ + S)\ 
where we have parametrized S as: 

{K + i)U + {K 



8iK U dU 
- 1 (C/2 - 1)2 



S = iK 



{K-i)U -{K + i) 



The variable U has been defined to send the poles S = —K"^ and 5 = — 1 to 
U = ±1. One gets the two parameters solution (parameters K and E) up to 
an origin for time, which we fix by requiring that i = for [/ = 0: 



t 



t-u 



or 



1 



tc 



c/2 



t-tc 



1 - C/2 ^ 



1 AiK 
~i^2 _ 1 



UJ 



We shall soon see that t = too is a second singularity of the dynamical variables, 
that we can express explicitly. For ease of comparison with the following, we 
present x±{t) = x{t) ± iy{t): 



x+ 



2Kg [{K-i)U-K-i] [{K + i)U + K-i] 1 

U 



(i^2 - l)2([/2 _ 1)2 



2Kg [{K -i)U -K-i] [{K + i)U + K 



(K2-l)2([/2_l)2 



2Kg [(K -i)U- K -i] [{K + i)U + K ■ 



A 



(K2 - 1)2(^72 _ 1)2 
(i^2-l)2(i^2^1)([;2_i)5 



U 



64K2 [{K -i)U -K-i] [{K + i)U + K-i] 
In terms of the t variable, we get the simpler expressions: 



x+{t) 



x^{t) 



2Kg 



^2(^2 _ 1)2 
W2(K2 _ 1)2 



{K^ + I] 



t-tc 



3/2 



1/2 



t - tr 



{K' + l)[ — 



3/2 



t - tr 



1/2 



4iK I — 



(26) 



(27) 



We see that x+ behaves as t 2 and x- behaves as when t ^ 0. If 

1 

we expand around t = we get Puiseux expansions in 1 2 . These expansions 
depend on three parameters, K and E plus the origin of time to- This is because 
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we are analyzing the trigonometric solution which fixes one of the constants 
to I = 0. We shall see later on that it can be generalized to a four parameter 
expansion in the elliptic case. The energy parameter appears factorized in front 
of x+ and x_ in the form oigjuP' = 2E/g. 

3 

Around ioo, we see that behaves as (t — ioo)^ and x_ behaves as (t — 
ioo)~^ which is symmetrical with the behaviour at t = 0. This is compatible 
with the fact that the equations of motion admit a symmetry {xj^{t)^x_{t)) ^ 
{-x-{t),-x+{t)). 

Remark that x±{t) are defined on the two sheeted covering of the Riemann 
sphere with two branch points at t = and t = too- The variable U that we 
have introduced is in fact a uniformizing variable for this covering, so that x±{t) 
are rational functions of U. Moreover U <-> —1/U corresponds to t ^ (too — t) 
and exchanges x+ and — x_. The extra minus sign means that we have to 
change the determination of the square root in the t variable. The U variable 
makes this completely unambiguous: 

=-^-(C^), z(^-^^=z{U), a(-1)=A(C/) 

We emphasize that, although the system is integrable, the solutions diverge 
with square root singularities at finite times t = 0, and t = too- 

We now return to the elliptic case. Let us define the variables X = S,"^ — 
E/{6g) and Y = rj^ - E/{6g)- The equations (1221231) become: 

1 f"" {X + E/{6g))dX , 1 {Y + E/{6g))dY 

t + CE = -T—zi I , , , h 



CI 

where now: 



Introducing the Weierstrass functions 

X = pi(Zi) = p(Zi,52(/),53(/)), Y = P2(Z2) ^ p(Z2,52(-/),ff3(-/)) 



1 dx 1 dY 



^i^J ^P+{X) ^i^/g J yfPJj) 



P±{X) = _ - 33 (±1) 

1 ^^2 , „ E / 2 9 



the above integrals reduce to: 

1 

t + CE ■- 



^{Zi + Z2) - Ci{Zi) - C2{Z2] 



(2S 



CI = ^[Z.-Z2] (29) 
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where ( is the Weierstrass zeta function, (' = —p. The p function has two 
periods 2u;j, j = 1, 2, so that p{z + 2ujj) = p{z), but the zeta function is quasi 
periodic, C,{z + 2ojj) = C,{z) + 2r]j. Here we have two set of periods ujj and 
r]j according to the function pi or p2, which are in fact functions iOj{±I) and 

Note that x±{t) have poles and zeroes when ^ ± irj vanish, that is when 
^'^ -\- r]"^ = X + Y + E/{3g) = 0. Hence we have to solve: 

E/{3g)+pi{Z^) + P2{Z2)=0 (30) 
Z1-Z2- 4i^ci = (31) 

But differentiating eqs. (|28|29p we find 5Z2 = bZ\ and 

<5t = ^ + pi(Zi) + P2(Z2)) '^^1 + ^ (P'l(^l) + P2(^2)) (5^l)' + - • • 



The first term vanishes when ^ ^t]^ = hence around such a zero bZ\ ~ \fbi. 
As a consequence, in view of eq. (fT8]) . x±(t) behaves as either or bi^l"^ at 

such a point, according to the vanishing of ^ + ir? or ^ — if]. Note this is similar 
to the trigonometric case. 

However finding the pattern of these singularities is messy, because in the 
equations (|30|31|1 we have two incommensurate lattices of periods for the two 
Weierstrass functions. However we can easily see that there is an infinite num- 
ber of singularities. This is because since the two lattices are incommensurate, 
for any large R and small e, one can choose Y in the first lattice and W 
in the second, such that |F — W\ < e and \V\, \W\ > R. Starting from a 
solution Zi, Z2 of our equations, we set Z[ = Zi + V and Z2 = Z2 + W , 
which still obey eq. ([30]) . However eq. ([3T]) is violated at order e. Choose 
Zl = Z[, Z2 = Z'l — M^ci and plug this in eq. ([30]) . It then gets of order e 
but this is an equation for the variable Z\ which has, by complex analyticity, 
an exact solution close to this approximate solution. Taking larger and larger 
values for R one gets an infinite number of solutions. Around each of these 
solutions we have Puiseux expansions in the variable bt^l"^ . 



4 Kowalevski analysis. 

If the swinging Atwood's machine is an algebraically integrable system the 
dynamical variables can be expressed algebraically in terms of a linear motion 
on some Abelian variety, in particular all variables and time can be complexified 
at will. We may expect that, for general initial conditions, the dynamical 
variables will blow out for some (in general complex) value of the time 



11 



t. Around this value the dynamical variables should have Laurent behavior, 
hence one expects to find Laurent solutions depending on parameters (initial 
conditions) if the phase space is of dimension A^. In practice one searchs for 
Laurent expansions at t = (one fixes Iq = 0) so an admissible Laurent solution 
should have A^ — 1 parameters, that is 3 parameters for the example at hand. 

The Puiseux solutions we have found in previous section have the following 
singularity: x and y blow up but z ^ 0, hence + — > 0. This means that 
the singular solutions are such that the mass m goes to the origin but rotating 
faster and faster. If we expand x and y in negative powers of t there must 
be large cancellations such that + ^ 0. It is much more convenient to 
factorize x^ + y^ and have the cancellation between the two factors. Reminding 
that: 

x± = X ± iy 

the equations of motion are 

mx+ = —img + Ax^ 

mx- = img + Ax_ 

Mz = -Mg - \z 

z^ = x+x_ (32) 

The value of A is a consequence of these equations: 

mM z^ — + y(y — z) 



A 



M + m z^ 

where y = —i{xj^ — x_)/2. Let us remark that this system of equations is 
invariant under (a;_|.,a;_) (— x_, — a;_|_), in particular y and A are invariant. 
The system is also invariant under a similarity transformation: 

1 



x±(t) ^ ii^x±(tlii), z(t) ^ m^z(Vm), \{t) ^ — A(i/^) 



We first analyze equations (|32]) at the leading order. We thus look for 
solutions of the form: 

x+ = ax€ ^ , x_ = hxfi H , 

so that eq. (f32ll requires 

z = ci t 2 + • • • , = a\b\ 

At lowest order we then have: 

^ ^ mM aibijp - q)HP+i-^ + Ag{y - z) 
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Clearly equations of motion (f32l) require that A behave as \/t^ for solutions 
blowing out as powers. At first sight there are two ways in which this can 
happen: when the first term in the numerator is dominant, or when the second 
term is dominant. We can always choose p < up to exchange of x+ and x_, 
hence p < Q since we want to have at least one dynamical variable diverging. 
On the other hand 2: ^ so g is positive, hence y — z = 0{t'P). The first term 
is dominant when q < 2, and for p / g one has indeed A ~ When q = 2 

both terms are of the same order and for g > 2 the second term is dominant, 
so that A = 0{t~'^) which is not allowed. Hence we have basically only two 
cases to consider, either p < 0, g < 2 in which the integrable case studied above 
belongs (p = —l/2,q = 3/2), or the case — 2 < p < 0,g = 2, which, as we will 
see, covers more general values of the mass ratio M/m. 



4.1 Integrable case. 

Since p < 0,q < 2 we have p + q — 2 < {p, q, ^^), and we can neglect the term 
g{y — z) at leading order in the expression of A. We find, for p ^ q: 

mM{p — q)^ 1 
~ 4(M + m) i2 + 

Similarly the equations of motion for x± give: 

p(p-l) = ^(^^(p-.)^ = .(.-l) (34) 

so that {p — q){p + q — ^) = hence, since p ^ q, p < q and we have p + q — 1 = 0. 
Since by positivity in eq.jM]), p and q cannot belong to [0,1] this implies, 
together with p > p + q — 2 = —1 that: 

-1 <p <0, 1 < q <2 

Using p + q = 1 the mass ratio takes the form: 

M = —Ampq = m[{p — q)^ — 1] = m[{2p — Vf' — 1] 

and the mass ratio M/m is thus in the interval ]0, 8[. 

The integrable case corresponds to M = 3m, and falls into this analysis 
with: 

1 3 

^ = -2' ^=2 

These exponents are exactly those we have found in the exact solution of the 
elliptic integrable case. There are no other values of p in ] — 1, 0[ compatible with 
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integer values of the mass ratio M/m which could, according to [4], correspond 
to seemingly integrable behaviour. We thus consider, in the following, the 
integrable case M/m = 3. 

As noted above the second conserved quantity is given in polar coordinates 
for m = 1, introducing for convenience H2 = I V2/8, by: 

-^H2 = r^e—{rcos{e/2)) + -(rsin0)(r cos(0/2)) 
2v2 dt 2 

which reads in cartesian coordinates as: 



H2 = (xy - yx) — {z^ - zy) + gx^^ z{z - y) 

^z{z - y) dt 

Taking the square to eliminate the square roots, we get: 

2 



1 



-{xij-yxf 



^ z{z - y) 
We can setup an expansion in powers of \/t 



^{z^ - zyU +2gx{xy-yx)^{z^-zy)+g^x^{z^-zy) 



X- 



z 
A 



t"2(ai + 02*2 H 

+ 62*^ + •••) 
t^2{di+d2t^ + •••) 
r'^ih + ht^ + ■■■) 



aibi = dj, h = ^ 



We already know that 

Inserting into the equations of motion, we find the recursive system: 



( \ 










Bs+i 


ds+i 




Ds+i 


V ^.+1 1 




\ Ls+i J 



IC{s) 



m 



(^-l)(^-3) 
4 



-bi 



u 







-ai 






2d^ 



+ h 



—ai \ 
-hi 
di 

/ 
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The square matrix in the left hand side is called the Kowalevski matrix, and 
the vector in the right hand side is given by 



s-1 

j=i 

s-1 

Bs+1 = ^lj+ibs-j+i + img6s,i 
i=i 

s-1 

^s+1 = - ^ lj+ids-j+1 - Mgdsfi 
i=i 

s-1 s-1 

Ls+i = — ^ dj+ids-j+i + ^ aj+ibs~j+i 
i=i i=i 

The determinant of the Kowalevski matrix reads 

2 j2 

det(/C(s)) = -^^(s + 2)s^{s - 2) 

It has a double zero at s = and a third zero at the integer value s = 2. Hence 
potentially three arbitrary constants may appear in the expansion. Indeed 
the miracle happens at the third level where the equations determining the 
coefficients 03,^3 are degenerate, leaving one extra constant 63 = ci. The 
rest of the expansion is then completely determined at all orders. We find in 
particular: 

_ dl idjg ScidjVi {Aicidf - 7bldi) gt 



{{2cidl + ibldi) c/2 + 12 61 c? ) t 



8bf 

, 3 iqt'^ 5 (2icidi- 6?) gt^ 



2 5 61 dl 

{{6cidi+3ibl) 52 _ mbicfdi) tl 
40 6f dl 



+ ••• 



The existence of such a "miracle" is exactly what S. Kowalevski noted in [8] 
for her integrable case of the top. For this to happen one needs that the 
determinant of K,{s) vanishes for the correct number of integer values of the 
recursive variable s, which allows for a new indeterminate to enter the expan- 
sion. Moreover in this case the linear system has to be solvable which is far 
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from guaranteed. The general solution of the equations of motion must ad- 
mit a power series expansion, which thus must depend on 2N — 1 arbitrary 
constants for a system of degrees of freedom. In our case we find a solu- 
tion depending correctly on three constants, which extends the trigonometric 
solution described above. 

Inserting these expansions into the formula for the energy (fTTl) we obtain: 

Similarly, the second conserved quantity reads: 

It is interesting to compare these general results to the expansion in the 
trigonometric case eqs. (|26|27l) . One finds: 

^ gjK^ + 1) 

g^VK'^ - l(K^ + 1) 

* 3 

32K2 

^.g^/K{K'^ + 1) 
cji(i^2 _ 1)1 

With these values one checks that = as it should be in the trigonometric 
case, and that H is indeed equal to E. 

The dynamical variables (x, y, z) and their time derivatives are expressed 
m power series of Vt. These power series have a non vanishing finite radius 
of convergence (we know this at least in the trigonometric case from the ex- 
act solution) and we can check it numerically. To do that we compute the 
d'Alembert quotient \an+i/an\ relative to a series On^" which tends to the 
inverse of the radius of convergence of this series when it exists. We present 
the result of this computation for high order n for the series 
and X{t) in the figure 

In this and following similar computations, all values are calculated with 
absolute precision rational numbers using a formal computation tool. This 
ensures accuracy of the result. 

Since the Kowalevski expansion converges in a disk, the parameters ci, di) 
appearing in these series, and the origin of time to, can be considered as co- 
ordinates on an open set of phase space near infinity [12]. The question then 
arises to compute the Poisson brackets in these coordinates. 



Cl = 

di = 
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d Alcmbcrt criteria. Intcgrablc case; cl = 1 



^\J\AA^ ' 



Figure 2: d'Alembert criterium for convergence for p=-l/2, q=3/2. 

To do that, we start from: 

{A,{t),x±{t)} = ±ix±{t) (35) 

This equation is vahd for any time since the time evolution is a canonical 
transformation. We thus insert into it the series for x±{t), where these series 
are really series in [t + to)2 . Similarly 

^^zit) = i^{x+x^ - x^x+){t) 

is expressed as a series in (t + to)^ and eq. ([35]) is an identity in t. The Poisson 
bracket is computed with the rule: 

ZdFdG dGdF\ , , , / dF dG dG dF\ , 

dF dG dGdF\ , f dF dG dG dF\ 

{to,di]+[——-——] {5i,ci} + 



dto ddi dto di J ' \ dbi dci dbi dci 

dF dG dG dF\ [ dF dG dG dF\ , 

{budi} + T^TTT - {ci,di} 



dbi ddi dbi ddi I \ dci ddi dci ddi 
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Plugging F = Az{t + to) and G = x±{t + Iq) and identifying term by term in 
(t + to) we get an infinite system for the six Poisson brackets of the coordinates, 
which is compatible, and whose solution is given by: 



{to,di} = 
{to,bi} = 



{to,ci} 



Amd\ 
bi 



{61, di} 



2mdi 

9^ + I661C1 
+ 326ici 



{ci,di} 



{ciM] 



32mdj 



We can then check that 



{H,h} = {H,ci} = {H,di} = 0, {H,to} = l 



Finally we see that canonical coordinates can be chosen to be the pair of cou- 
ples {H,to) and {logbi, mdl), hence the Kowalevski constants are essentially 
Darboux coordinates in a neighbourhood of infinity. 

This shows the interest of these Darboux coordinates in a vicinity of infinity, 
but the whole question of integr ability is a global one. Our problem is therefore 
to try to extract some information from the Kowalevski series beyond their 
disk of convergence. In the following we investigate this problem numerically. 
First we have seen that o„+i/o„ tends to a complex number that we call with 
hindsight to^^"^ . Hence a„ behaves asymptotically as a„ ~ too"^^- One can do 
even better and look at the prefactor. Assuming that 

we can extract the coefficient a by computing the quantity: 



We show the result of this calculation in figure jS]). Note that the curves begin 
by large oscillations but for n sufficiently large, in the asymptotic regime, the 
exponents a tend to constants. Comparing with the dominant terms in the 
binomial formula: 



lim n 



2 0,n—20-n 



— 1 = —a 






we see that setting z = \J t/t, 



we read from figure ([3]) the various exponents: 
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x+{t) 


~ (1 - 




x.{t) 


~ (1 - 




z{t) 


~ (1 - 




m 


~ (1 - 





The consequence of this observation is that x±{t) have Kowalevski expan- 
sions around too with indices which are exchanged as compared to those around 
t = 0. Hence we know that: 

x+ = -b[ (t^-t)!- '^ ~ _ (too - t)i + • • • 

where we have introduced a change of sign required by the symmetry x± 
— x^, and the symmetry of the equations of motion under t ^ too — t. The 
series expansions have new parameters b'^ , c'l and d'l . In the trigonometric case 
we see from the exphcit formulae that they are equal to the original parameters, 
see eqs. (l26l[27ll . 
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We have learned from the previous analysis that the singularities are al- 
ways of the Kowalevski type, with well-defined exponents. This is perfectly 
consistent with the exact solution in the trigonometric and elliptic case. 



4.2 Non integrable case. 

We now explore the region of parameters —2<p<0,q = 2. We assume that: 



x+ ~ ait^, x_ ~ z ~ cita"*"^, cf = aib 



Notice that z — > since we assume p > —2, and that y = — — x_) ~ 
— ^ait^. We see that both terms in eq. (f33]l for A contribute: 



^ mM ((P_.Y_ i9 \ ]_ 
^ M + m[\2 J 26Jt2 



The x± equation give: 



, mM f fp \2 ig 

2mbi = imq H 1 — 

^ M + m\\2 J 2bi 



Solving for bi we find 



M = -4m^, b, 



p + 2' {p-2){p + l) 

Notice that the mass ratio is positive if — 2 < p < 0, and that: 

mp{p — 1) 
A ~ 

For relatively prime integers r and k we set: 

r 

p = — — , —2k < —r < —k 
k 

We perform the Puiseux expansions: 

x+ = t (ai + 02** H ) 

X- = t^{bi + b2t^ +■■■) 

Z = + (i2t^ H ) 

A = t-^ih + ht^ + ■■■) 
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We already know that 



r{r + k) dl 



^ , M = 4m- 



2k -r 

When we plug this into the equations of motion, we get a system of the form: 



where the Kowalevski matrix reads: 



( ^s+l \ 










Bs+1 


ds+i 




Ds+i 


\ls+lj 




\Ls+iJ 



(36) 








^ - h 








j^^ ik+s-r/2)is-r/2) 

\ —bi —ai 2di 

and the right hand side of equation Eg is given by: 



-ai 


\ 


-hi 




dl 







/ 



Ag+i = ajk+2-j - img6s,2k+r 

3=2 



B 



D 



L 



s+l - X] ^3^s+2-j 

s 

s+l = - XI djls+2-j - Mg6s^k+r/2 

s s 

s+l = — X djds+2-j + X ajbs+2- 

j=2 j=2 



(37) 



For s = 1, the quantities A2, B2, D2, L2 are meant to be zero. The determinant 
of the Kowalevski matrix reads: 

det(/C(s)) = -6m^dlj^^^s{s + k){s - r){s + k - r) 

In order that this determinant vanishes for two positive integer values of s, 
assuming A; > 0, we should have 



r > 0, r -k> 



21 



Prom the equation for D^+i it is natural to choose r even, otherwise the weight 
Mg would disappear from the problem which is not physical. In order for r/k 
to be irreducible, we must choose k odd. Setting r = 2r' we finally get: 

|<r'<fc, P = -2j (38) 

When things are setup this way the Kowalevski determinant has two strictly 
positive integer roots, so that, potentially three arbitrary constants enter the 
expansion, or the expansion is impossible. Impossibility occurs when the right- 
hand side of equation Es is non vanishing and doesn't belong to the image of 
JC{s) for values of s which are Kowalevski indices. It turns out that in most 
cases the right-hand side vanishes as we now show. 

First, since we want to examine the behavior for s = r — k and s = r we 
can limit ourselves to studying the system for s = 1, • • • , r. In this case the 
Kronecker deltas in eqs. ([371l always vanish. For 5s,2fc+r it is obvious, for Ss^f^^j./2 
note that, since k > l + r/2 we have k + r/2 > r + 1. Since the induction 
starts with A2 = B2 = D2 = L2 = 0, we get, if s = 1 is not a Kowalevski 
index that 02 = 62 = ^2 = ^2 = 0, hence the right-hand side vanishes for 
the next equation s = 2. This goes all the way up to s = r — fc, hence when 
we hit the first Kowalevski index, it is always with vanishing right-hand side. 
The existence of a non trivial solution a^-fc+i, • • • , Ir-k+i is thus guaranteed. 
Let us assume for the time being that the first Kowalevski index is such that 
{r -k) > 1, that is r > A; + 2. 

As a consequence of this previous step, when s = r — A: + lwe find that 
As+i reduces to aj,+fc_i/2 which also vanishes because I2 = 0. More generally 
we have As+i = Yli]=r-k+i %'^s+2-j which vanishes when s + 2 — j < r — k + 1 
for all j in the sum, and similarly for the other components. This occurs when 
s < 2{r — k). For s = 2(r — k) the right-hand side of equation Eg doesn't 
vanish, and assuming we are not on a Kowalevski index, there is a unique non 
vanishing solution a2{r-k)+ii ' " , h{r-k)+i- The process continues and it is easy 
to show by induction that the right-hand side of equation Eg doesn't vanish 
only for s = n{r — k),n positive integer, so non trivial solutions are of the form 
^n{r-fc)+i) ■ ■ ■ ) ^n(r-fc)+i- Indeed, to get a non vanishing ajls+2-j we need to 
have j = n(r—k) + l and s+2— j = n'{r—k) + l so that s = {n+n')[r—k). In this 
case only we have A^+i, • • • , -Ls+i and thus a^+i, • • • , Ig+i, non vanishing. This 
shows that the next non vanishing positions are of the form (n + n'){r — k) + l 
establishing the recurrence. 

The second Kowalevski index is s = r and this cannot be of the form 
n(r — k). Indeed, since r and k are relatively prime, if we have r = n(r — k) 
we get (n — l)r = nk, hence n = pr and n — 1 = qk for some integers p and 
q. Then (n — l)r = nk = qkr = prk so that p = q and finally 1 = p{r — k) 



22 



which is only possible for p = 1 and r = k + 1. This is precisely the case we 
have excluded up to now. As a consequence, when we arrive at the second 
Kowalevski index s = r, the right-hand side of equation Sg vanishes and there 
is a non trivial solution, with an extra constant. 

We have shown that two new constants of motion always appear for all cases 
r = k + 3,k + 5, ■ ■ ■ ,r = 2k — 2. This covers an infinite number of values of 
the mass ratio M/m for which the Kowalevski criterion is satisfied (with weak 
Painleve solutions), but for which the system is presumably non integrable. 

Finally we discuss the case k = r+1. The first Kowalevski index is s = 1. In 
this case the right-hand side vanishes and we have automatically a non trivial 
solution [^2, ^2, ^2; '2]- From this point, all other solutions of the linear system 
don't vanish, and in particular, for the second Kowalevski index, s = r, the 
right-hand side of the system is not trivial. For a solution to exist it must be in 
the image of /C(r). Equivalently, let us consider a covector U = [ui,U2,U3,U4] 
such that C/./C(r) = 0. Explicitly: 

ui = 2g'^k^, U2 = dl{k + l){3k + lf, U3 = i di g k"^ (k - l){3k + 1) 

ua = -2imgk{k + l){2k + l){3k + 1) 
The condition to be satisfied is that the scalar product: 

W{s) = UiAs+i + U2-Bs+1 + U2,Ds+l + UiLs+l 



of this covector and the right-hand side of eq. (l36l) vanishes for s = r = A; + 1. 

For arbitrary k and s = 3, 4, ... we have computed this scalar product W{s), 
and we have observed that W{s) has a factor {s — k + 1). For example we get: 

W{s = 3) = -I^M^ - + 1)(2^ + + 1)' 



VA:5(A; + 2) 

Note the factor {k — 2) = [r — s). For s = 4 we next get: 

imcfdjik - 3){k + l){2k + 1)(3A; + l)^P6(fe) 



W{s = 4) 



9653(A;- l)2fc8(fc + 2)2(A; + 3) 



with the factor {k — 3) = (r — s). Here ci is the Kowalevski constant which has 
be introduced at s = 1, and Peik) is some polynomial in k of degree 6. The 
factors in the denominator of course come from similar factors in det(/C(s)). 
The expression for s = 5 has the same type of factors in the numerator and 
denominator, with a more complicated polynomial Pjik) and always a factor 
(r — s). This behavior is persistent as far as one can compute. The consequence 
of the presence of the factor (r — s) is that, for any k, when we arrive at the 
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second Kowalevski index, s = r = k + 1, the scalar product W{k + 1) vanishes 
and the Hnear system is solvable. We can thus state that for all admissible pairs 
{k, r) the swinging Atwood machine has weak Painleve expansions depending 
on the full set of parameters. 

For example an interesting case occurs when the mass ratio M/m = 15 
where the system doesn't look chaotic, see [2]. This case is obtained when 
k = 19 and r = 26. The linear system is solvable in this case, although the 
new arbitrary constants occur very far from the beginning of the expansion. 
We shall refrain to exhibit the solution in this case, since it is very bulky, and 
proceed to show what happens with smaller values of k and r. 



4.3 Example: the case k = 3, r = 4. 

When k = 3 we have necessarily r = 4. The Kowalevski exponents are s = 0, 
s = 1, s = 4. The dynamical variables x± expand in Puiseux series of t^^^ 
which take the form: 



X4 



4 ,9 fWi 1 140«c? 2 14000c? 3 
^\9g 729 59049 



1960 z cf m- 32805 ic2r * 

H 5 —t3 H 

91854 c/5m 

n / 9 i o 1 7 i c? 2 14 c? 3 

X- = +cit3 + —^ts + —l-ts + 



10 30 g 243 g^ 

(96124 icfm- 177147 i C2 g'^) ^4 ^ 



918540 g^ 



m 



This solution depends on 4 arbitrary constants: to,di,ci,C2 (in the above ex- 
pansions t should always be understood as t + to)- We obtain 

- _ 5 ^1 (13412 cfm- 19683 C2 g^) 

91854 5^ ^ ^ 

The above constants can be used as local coordinates on phase space. To 
compute the Poisson brackets of the Kowalevski constants, we proceed as in 
the previous section considering {Az{t),x±{t)} = ±ix±{t). We find: 

{to,ci} = (40) 

{to,di} = (41) 
14 1 

{fo.e.) = (42) 
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D Alcmbcrt criteria, Non intcgrabic case: cl= 1 , c2= 2 



100 150 200 250 300 350 400 450 



Figure 4: d'Alembert criterium for convergence, liman+i/dn in the non integrable 
case ci = 1, C2 = 2 for = 450. 



1 • 13412 cf 



. (13412 cfm- 19683 C2 5^) . 7^ 
^^^'^^^ = 65610 dlg-^m = ''25^4 ^'^^ 

It is remarquable that these six relations ensure the compatibihty of an infinite 
set of relations. One verifies easily the Jacobi identity in spite of the crazy 
numbers appearing. We can compute the Poisson brackets with H 

{HM = 1 
{H,di} = 
{H,ci} = 
{H,C2} = 

so that to is the conjugate variable of H as it should be and the other ones 
are constants of motion. Notice that ((if,ci) is a pair of canonical variables 
commuting with the pair {H, to)- Kowalevski constants are essentially Darboux 
coordinates. 
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If there were an extra conserved quantity it would therefore be a function 
F(ci, C2,di). The variable C2 can be eliminated through H so that we can write 
as well F{H, ci,di). 

As in the integrable case we can compute numerically the radius of conver- 
gence, and the exponents which nicely fit with the above Kowalevski analysis, 
as shown in Figure ([5]). 

Exponents. Non integrable case: c1= ! , c2= 2 



x+ 

z 





100 150 200 250 300 350 400 450 

n 



Figure 5: Exponents at singularities in the non integrable case Ci = 1, C2 = 2 
= 450. 



To go further, we also compute the Fade approximants of the series. It 
is more convenient to consider the logarithmic derivatives x±/x± because the 
residues of the poles are the exponents. We present the polar decomposition 
of the [74,75] Fade approximant fo This shows clearly eight true 

singularities with residues respectively -1.33 and 2 (up to numerical errors) 
consistent with the Kowalevski analysis. The other poles having small residues 
correspond to strings of poles and zeroes representing algebraic branch cuts in 
the Fade analysis. Note we have set t = and we have cancelled the leading 
z"^ at the origin. 

2.07 .0366 i .0295 0.016 i 

^+'^+ ~ 0.812 + 0.618 f -Fz ^ 0.813 + 0.622 i -Fz ~^ 
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Log Padc, Not! Intcgrablc case. cl = 1 c2= 2 N= 1 19 
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Figure 6: Poles and zeroes of Pade approximant [M,M + 1] of x^/x^ in the non 
integrable case ci = 1, C2 = 2, M = 59, and = 119. 



2.05 - .0136i .0351 - .00792 i 

"^-0.33 + 1.04i + z ^ -0.332 + 1.04i + z " ' 
2. 13 -. 0627 i .0725 -. 0176 f 

"^-0.95-0.012^ + 2 ^ -0.954 - 1.33 x lO^^^ + z 
-1.34- .00154i -.00711 + 0.0025 i 

"^-0.637- 0.83 X i + z -6.49 x 10"! - 8.52 x 10"! i + z 
2.38- 0.027i .0281 + .0125i 

"^-0.192 - 0.703 i + z ^ -0.192 - 0.705 i + z " ' 
-1.34 + 3.547 X 10"^ i -0.0064 - 3.344 x 10"^ i 

0.175 - 0.84i + z 0.177 -0.85i + z 

2.29 + . 114 i 0.077 -. 0335 i 

"^0.629 -0.545i + z 0.63 - 0.547 i + z 

-1.34- .00401 i .0802 -. 0866 i 

"^1.45 - 0.586 X W-^i + z ^ 1.62 - 0.77f + z "^ " ' 

We see that this structure is very similar to the one we have observed in the 
integrable elliptic case. This semi-local analysis doesn't appear to be able to 
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discriminate between the integrable and non integrable cases. 

5 Conclusion. 

We have studied the swinging Atwood machine, which is beheved to be non 
integrable except for the mass ratio M/m = 3. We have shown on the expHcit 
solution of the integrable case that the Kowalevski analysis is valid, but requires 
weak Painleve expansions. We have extended this weak Painleve analysis for 
other values of the mass ratio, and shown that it is valid for an infinite number 
of cases. Hence this model is remarkable in that it exhibits an infinite number 
of cases where the Kowalevski analysis works at the price of using Puiseux 
expansions. However only one of these cases is known to be integrable, while 
the other ones are believed to be not integrable. 

In the cases where Kowalevski expansions are available, we have shown that 
the constants appearing in these expansions provide Darboux coordinates on 
an open set of phase space around infinity. The question of integrability of 
the system therefore reduces to the global nature of this coordinate system 
(to, Cl, C2, di) on phase space. 

On this open set, knowing the Poisson brackets eqs. (|40ll45]l . we can try to 
find the conjugate variable of Iq. We find that H must be of the form: 

15 

H = -Y^djc2 + h{ci,di) 

The first term agrees with the exact formula in equation ([391) . The function 
h{ci,di) is not determined but it is of course crucial to have a "good" function 
H{x+,X-,x+,X-). Clearly we can, in principle, invert locally the system of 
equations 

X+ = X+{t -to,Ci,C2,di) 
X- = X-{t -to,Ci,C2,di) 
X+ = X+{t - to,Ci,C2,di) 
X- = X-{t - to,Ci,C2,di) 

where in the right hand sides we mean the Kowalevski series. In doing so, we 
will find 

t — to = r(x+, x_, x_) 

Cl = Cl{x^,X^,X-^-,X^) 
C2 = C2{x+,X-,X+,X-) 

di = Di{x+,X-,x+,x^) 
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but the functions T, Ci, C2, Di will behave in general extremely badly. All this 
shows that it is in general impossible to make statements about the integrability 
of the system on the only basis of the Kowalevski analysis. In this context it 
is remarkable that the global hamiltonian indeed exists, and it is even more 
remarkable that a second global hamiltonian exists in the integrable case. We 
see here in a striking way the global nature of integrability. 

In the non integrable case, in an attempt to progress beyond the analysis of 
a single singularity, we have used Pade expansions. In this semi-local analysis, 
the panorama which appears is still remarkably similar to the one appearing 
in the elliptic integrable case. Hence Kowalevski analysis is not sufHcient to 
characterize integrability. Nevertheless it is a very non trivial property whose 
significance remains mysterious. 
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